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Цель работы заключается в повышении быстродействия 
устройства обращения ковариационной матрицы помех 
адаптивной антенной решетки за счет сокращения числа 
выполняемых операций. Это достигается использованием на 
этапе разработки алгоритма обращения априорной инфор- 
мации о свойстве эрмитовости обращаемой матрицы. В 
отличие от известных алгоритмов обращения, базирующих- 
ся на применении метода Гаусса — Жордана, в основу 
предложенного алгоритма положен метод окаймления. 
Актуальность разработки обусловлена сложностью метода 
Гаусса — Жордана и необходимостью большого числа 
операций при его использовании. Указанные особенности не 
позволяют реализовать режим реального времени при обра- 
ботке сигналов в вычислительных устройствах адаптивных 
антенных решеток, широко применяемых в системах связи, 
радиолокации и радионавигации. Предложенный метод, 
дополняющий известный метод окаймления учетом свойств 
эрмитовости ковариационной матрицы помех, позволяет 
построить алгоритм на базе рекуррентных соотношений. 
Получаемый при этом выигрыш от сокращения объема 
вычислений составляет не менее 25 % по сравнению с мето- 
дом Гаусса — Жордана. Уменьшение объема вычислитель- 
ных затрат, а также более простой вид соотношений, приме- 
няемых для построения алгоритма обращения матрицы, 
дали возможность разработать и более простую схему 
устройства, которое можно использовать в процессорах 
адаптивных антенных решеток для получения обратной 
матрицы. 


Ключевые слова: адаптивная антенная решетка, вычисли- 
тельный блок адаптивной антенной решетки, обращение 
ковариационной матрицы помех, метод окаймления, свой- 
ство эрмитовости ковариационной матрицы, сокращение 
объема вычислений, устройство для реализации процесса 
обращения матрицы. 
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Введение. На современном этапе развития радиоэлектронных систем (РЭС) в области связи, радиолокации и радиона- 
вигации отмечается значительное усложнение электромагнитной обстановки. Это связано с высокой пространствен- 
ной плотностью размещения РЭС и ограничениями используемых частотных диапазонов. Так, по данным [1], число 
базовых станций формата 3С/4С только одного российского оператора «МегаФон» к концу 2013 года составляло по- 
рядка 43,5 тыс. 

В данных условиях, как правило, параметры мешающих сигналов (в первую очередь, направление прихода), 
неизвестны. Поэтому для обеспечения помехоустойчивого приема применяются методы адаптивного формирования 
«нулей» диаграммы направленности (ДН) в многоканальной антенной системе — адаптивной антенной решетке 
(ААР) [2-4]. Функциональная схема получаемой при этом адаптивной антенны показана на рис. 1 [2]. 
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Рис. 1. Функциональная схема антенны 


Конструкция антенны (схема расположения излучающих элементов) выбирается в зависимости от требуемых 
характеристик. В общем случае схема включает многоэлементную антенную решетку, диаграммообразующую схему 
и цифровой адаптивный процессор. Антенная решетка используется для приема сигналов, которые суммируются с 
весовыми коэффициентами, вычисленными_адаптивным процессором в диаграммообразующей схеме для формирова- 
ния диаграммы направленности антенны. 

Обработка сигналов в адаптивном процессоре базируется на использовании современных методов спектраль- 
ного анализа [2-4]. Они позволяют определять число и угловые координаты источников излучения, не прибегая к 
электрическому или механическому перемещению диаграммы направленности антенны. Используются исключитель- 
но алгоритмические способы обработки сигналов, принятых элементами антенной решетки. В результате можно в 
режиме реального времени отслеживать координаты источников излучения, находящихся в зоне наблюдения. В этом 
заключается основное достоинство методов спектрального анализа. Общая их особенность — использование взаимно- 
спектральной матрицы сигнала (ковариационной матрицы), оцениваемой на некоторой заданной частоте или в неко- 
торой узкой полосе частот выходного сигнала приемника. С точки зрения теории, методы локализации источников 
излучения могут быть применимы только с некоторой задержкой, связанной с накоплением определенного числа 
мгновенных выборок (снимков) выходов приемников. Причем чем больше размер накопленных выборок, тем точнее 
получаются результаты. 

Используемые в адаптивном процессоре алгоритмы базируются на нелинейной обработке сигналов методами 
линейной алгебры и теории матрицы с обращением или спектральным разложением корреляционной матрицы вход- 
ных сигналов, имеющей большую размерность. Реализовать такую обработку в реальном масштабе времени в насто- 
ящее время затруднительно даже на самой современной элементной базе [5]. В связи с этим актуальной является зада- 
ча поиска алгоритмов, упрощающих процесс обращения с учетом свойств матрицы, включая и алгоритмы получения 
решения в аналитической или полуаналитической форме [6, 7]. Так, в ряде работ, например [8-16], рассматриваются 
алгоритмы формирования обратной матрицы. В работе [17] описывается устройство, реализующее процесс получения 
обратной матрицы. Его построение основывается на использовании метода Гаусса — Жордана, представленного ре- 
куррентными соотношениями. В работах [15, 16, 18] предложен алгоритм обращения, основанный на использовании 
метода окаймления [19]. Однако в указанных алгоритмах не учитываются свойства обращаемой матрицы. 

Цель работы — повышение быстродействия устройства обращения ковариационной матрицы помех адаптив- 
ной антенной решетки путем сокращения числа выполняемых операций на основе использования в алгоритме обра- 
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щения априорной информации о свойствах эрмитовости матрицы. 
Выявление свойств ковариационной матрицы. Для выявления свойств обращаемой ковариационной матрицы по- 
мех рассмотрим модель задачи определения направления прихода сигнала. 

Математическое описание алгоритма определения направления прихода сигнала формулируется следующим 
образом [2-4, 8]. Пусть М№-элементная антенная решетка с идеальной точностью принимает сигналы, передаваемые 
некоррелированными источниками численностью М. При этом пусть передающая среда является изотропной и не 
рассеивающей, а источники излучения находятся в ее дальней зоне. Пусть канал системы передачи является узкопо- 
лосным. Значит, по мере прохождения с направления 9, (т= 1,2... М) радиосигнала $(0 =[5, (0,5, (0,....5и (0 
через решетку в виде волновых фронтов его огибающая остается неизменной. Здесь Т — символ транспонирования. 
Предположим, что отклик решетки на каждый сигнал является функцией только одного углового параметра 0,„. Тре- 


буется определить направление прихода волны 6,. 

С учетом сделанных предположений пространственно-временная обработка сигнала разделяется на простран- 
ственную и временную, выполняемые в произвольном порядке. В результате суммарный измеряемый сигнал опреде- 
ляется выражением [2]: 

х(1) = 4(0)5(+п(0, (0 
где и(®) =[и (@), идете — комплексный низкочастотный сигнал, эквивалентный полученному сигнальному 


вектору на антенной решетке в момент времени #; 4(0) =[а(@},а(0, ),..-- (9 М Г — вектор отклика решетки в направле- 


нии 0,. 

Известно большое число алгоритмов для обработки принятых антенными элементами данных с целью опре- 
деления направления на источник излучения [2]. Так, для алгоритма «сверхразрешения», подробно описанного в [2, 
4], направления прихода сигналов определяются положением пиков в пространственном распределении мощности: 


Р=(А"(@)КыА(6)) (2) 
Здесь корреляционная матрица входных сигналов определяется выражением 
В. = А(@)В.А"(@) + о7Г, (3) 
а функция ковариации входных сигналов определяется соотношением 
В=В ото эх (Е) (4) 
2 М 247» " 


где № — размер выборки; о” — дисперсия шума; / — единичная матрица; Н — символ эрмитового сопряжения. 
Анализ соотношений (2)-(4) показывает, что для получения информации о направлениях прихода сигнала 

необходимо осуществить обращение ковариационной матрицы. Ее размерность а также точность получаемых резуль- 

татов определяются размером выборки: чем больше объем выборки, тем более точные результаты могут быть получе- 

ны. Обращаемая матрица, как несложно заметить из анализа выражения (4), обладает свойством эрмитовости [19]. 
Таким образом, при разработке алгоритма обращения корреляционной матрицы помех, имеющей большую 

размерность, необходимо учитывать ее эрмитовость. 

Алгоритм обращения матрицы с учетом эрмитовости. Пусть эрмитовая матрица М порядка М№ является заданной 


и определяется следующим выражением: 


бп 512 тм 
5 5 ... 9 
21 22 2мМ 
М = (5) 
Ям 5м2 -* ММ 


В соответствии с предлагаемым алгоритмом на первом шаге из элементов 5,, (1=1,...,№М, Л=1,..,М№М) матри- 


цы М выбираются элементы 91, 50, 51, 5. Данные элементы позволяют сформировать блочную матрицу М." 


размерности 2х2. Формулы, на основе которых вычисляются данные элементы, — т), то, т, т® с учетом эрми- 
товости матрицы М имеют вид: 


(0 в о Ра 
Ти =| 511522 [52| 5525 ИП = [511522 [52| 512 


(__ — 2\-! = № 2-1 * _ (1. 
Ты = | 511522 [52] 51 =—| 511522 о 52 => › 
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На втором шаге происходит формирование обратной матрицы третьего порядка М.'. Для построения данной 


2\-1 
(лы вы] ) бт: 


матрицы дополнительно используются элементы $55, 553 533 5315) Ковариационной матрицы М, окаймляющие полу- 
ченную на первом шаге матрицу М.', как приведено ниже: 


(0 (1) 


ти тТЬ 5 91 м 
(1) (1) 
т ТЬ 53 52 
51 52 533 53мм (7) 
5м 5м№2 5мз 5 му 


© * * 
Однако для эрмитовои матрицы 5. =5.; , 5. =5.› . Поэтому для нахождения матрицы третьего порядка до- 
статочно использовать только элементы 5.1, 5.5, 533. Элементы этой матрицы находятся с помощью формул: 


(2) (2) (а) („а („а 


ти ть ти ть ти т я. ти т 
(2) = КО КО (55) 7 (55) Фа |? — 
т т» т т т >> 21 >> 
(2) (0 „а 
тз ту т, =т 1 9 
| | пФ ю (5315) И’ (9) 
т т» ть 
(0 „(0 
1 т’ т 
2 ты 
(тб т: = Я (55) пе (10) 
21 > 
1 
(2) _ 
тие, (0 
0 „0 
а ти т *т 
где = 53 — (51 52 ) [Е ь 


то 


(0 
21 т 


22 


—1 
п+1 


В общем случае для перехода от обратной матрицы М, | порядка и к обратной матрице М... порядка и+1, 


выполняемого на и-м шаге преобразования, дополнительно используются приведенные на рис. 2: 


5 


— элементы а а ‚ образующие матрипцу-столбец В; 


5 


2би-и 


— элементы вы } ‚ обозначенные как матрица-строка С’; 


— элемент 5 представляющий блок О. 


п-+1и-1 2 


д 


5 им 


5 7+1 № 








$ мм ) 





Рис. 2. Формирование структуры обращаемой матрицы 


С учетом эрмитовых свойств матрица-столбец В равна комплексно сопряженной и транспонированной мат- 


—1 
п-+1 > 


рице-строке С ‚то есть В =С"Т. Следовательно, матрица М... ‚ элементы которой определяются формулами, описан- 


ными в [19], примет вид: 
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и [Моем утеми —нтемтут 
М 1 = п п п п 


› 12 
п+1 -Н СМ” НН" ( ) 


где Н=р-СМ7С". 

Оценка эффективности предложенного алгоритма. Для сравнения эффективности предложенного алгоритма обра- 
щения с известным алгоритмом на основе метода Гаусса — Жордана воспользуемся результатами работы [20]. Анализ 
соотношений, описывающих предложенный алгоритм обращения матрицы, показывает следующее. Для обращения 
матрицы В „ размерности №х № с использованием метода Гаусса — Жордана потребуется (№)° операций умноже- 
ния и сложения. При использовании предложенного метода окаймления с учетом свойств эрмитовости корреляцион- 


ной матрицы В „ необходимо 0,75(№)’ +(№)? операций. Таким образом, предложенный алгоритм обеспечивает для 
матриц больших размерностей, когда № —> оо ‚ сокращение числа операций до (1-—(0,75-1/ №))-100 % = 25 %. 


Техническая реализация алгоритма. Схема устройства, реализующего предложенный алгоритм, отличается от 
устройства, описанного в работе [18] и основанного на использовании метода окаймления. При учете свойства эрми- 
товости обращаемой матрицы число вычислительных модулей может быть сокращено с 11 до 9, упрощается схемное 
выполнение вычислительных модулей, изменяются связи между ними. Кроме того, изготовление устройства обраще- 
ния матрицы для реализации предложенного алгоритма не представляет особых затруднений, поскольку в операциях 
матричного умножения используются типовые цифровые блоки, выполняющие операции умножения, сложения и вы- 
читания. 

Блоки хранения коэффициентов матрицы могут быть выполнены на основе адресного запоминающего 
устройства, включающего блок памяти, дешифраторы и формирователи с регистром адреса, формирователи записи и 
усилители считывания с подключенным к ним регистром числа, а также блок синхронизации (управления). 

Таким образом, устройство обращения ковариационной матрицы помеховых сигналов с учетом ее эрмитовых 

свойств позволяет за счет уменьшения числа вычислительных модулей упростить схему по сравнению с предложен- 
ной в [17]. Кроме того, с учетом сокращения требуемых для обращения операций можно повысить быстродействие 
работы устройства. 
Заключение. В цифровых адаптивных процессорах радиоэлектронных средств при определении направления прихода 
сигнала широко применяются алгоритмы обращения корреляционной матрицы помеховых сигналов на основе учета 
ее свойств. Это определяет актуальность разработки подобных алгоритмов. При большом количестве сигналов, обу- 
словленном плотным размещением средств радиосвязи в городских условиях, известные алгоритмы не обеспечивают 
вычисление пеленга в реальном режиме времени. Отличием предложенного алгоритма обращения является учет эр- 
митовых свойств ковариационной матрицы помеховых сигналов, позволивший почти на 25 % сократить объем вычис- 
лений. В результате устройство, реализующее операцию обращения корреляционной матрицы помех в цифровом вы- 
числительном процессоре, будет более компактным по сравнению с известными, а также обеспечит более высокую 
скорость вычисления пеленга сигнала. 
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